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S.  B.  Batdorf  and  H.  L.  Heinisch,  Jr. 

Abstract 

A primary  objective  of  statistical  fracture  theory  is  to  predict  the  prob- 
ability of  failure  for  an  arbitrary  stress  state  when  the  failure  statistics 
are  known  for  a particular  stress  state,  e.g.,  simple  tension.  It  is  proved 
that  the  rule  for  accomplishing  this  given  without  proof  by  Weibull  in  1939  is 
valid  for  shear- insensitive  cracks,  l.e.,  on  the  assumption  that  only  the  com- 
ponent of  stress  normal  to  a crack  plane  contributes  to  its  fracture.  Four 
different  failure  criteria  for  shear-sensitive  cracks  are  considered,  and  the 
results  are  compared  with  test  data. 
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Recently  Batdorf  and  Crose  proposed  a statistical  theory  in  which  atten- 
tion is  focussed  on  the  cracks,  and  their  failure  under  stress.  For  sim- 
plicity the  fracture  of  a crack  was  assumed  to  depend  only  on  the  component 
of  stress  normal  to  the  crack  plane.  The  theory  leads  in  a very  straight- 
forward fashion  to  a technique  for  finding  the  fracture  statistics  for  poly- 
axial  stress  states  when  they  are  known  for  the  uniaxial  case. 

A number  of  investigators  have  questioned  the  legitimacy  of  Weibull's 
rule  for  polyaxial  stress  states,  and  some  have  advocated  the  so-called 
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"postulate  of  independence".  ’ ’ Batdorf  showed  that  use  of  this  postulate 
is  generally  unconservative,  and  compared  its  predictions  with  those  of 
Weibull's  theory  for  biaxial  tension. 

It  is  shown  in  the  next  section  that  Weibull's  rule  for  treating  poly- 
axial  stress  states  is  equivalent  to  the  technique  derived  in  Reference  4. 

This  means  that  Weibull's  rule  is  correct  if  cracks  obey  the  assumptions  of 
Reference  4,  e.g.,  independence,  random  orientation,  and  shear  insensitivity. 
Since  the  presence  of  shear  on  the  crack  plane  reduces  the  normal  stress  needed 
to  produce  fracture,  it  brings  the  curves  for  uniaxial  tension  and  biaxial 
tension  closer  together.  Thus  if  uniaxial  data  are  used  as  the  basis,  biaxial 
fracture  stress  predictions  will  be  conservative  when  shear  on  the  crack  plane 
is  neglected.  However,  better  results  can  be  obtained  by  using  a better  frac- 
ture criterion. 

When  shear  is  taken  into  account  a complication  arises  in  that  it  is 
necessary  to  specify  the  shape  of  the  cracks,  and  sometimes  Poisson's  ratio 
as  well.  An  additional  complication  is  that  there  is  not  as  yet  a concensus 
regarding  how  to  treat  mixed  mode  fracture.  Because  of  this,  in  the  present 
paper  several  criteria  are  treated.  The  results  obtained  are  compared  with 
test  data.  Agreement  is  best  for  the  relatively  shear  sensitive  crack  models 
arising  from  strain  energy  release  rate  considerations. 


We  first  derive  the  fundamental  equation  of  weakest  link  theory.  Let  it 
be  assumed  that  a stressed  solid  can  fail  due  to  any  of  a number  of  independent 
and  mutually  exclusive  mechanisms  or  causes,  each  having  an  infinitesimal 
probability  of  failure  j . The  probability  that  the  i'th  mechanism  will 


■ — m 


not  cause  failure  is  ^sj  = 1 - ^APf j . The  overall  probability  of  survival 
is  the  product  of  the  individual  probabilities  of  survival,  i.e.. 


The  sum  of  the  individual  probabilities  of  failure  appearing  in  the  final 
equality  above  was  called  by  Weibull  the  "risk  of  failure"  and  was  given  the 
symbol  B. 

The  potential  causes  of  failure  are  the  individual  cracks.  For  purposes 

of  analysis  it  is  convenient  to  group  the  cracks  according  to  location,  the 

applied  stress  state,  and  crack  critical  stress.  We  assume  that  the  stress 

state  varies  slowly  so  that  within  a volume  element  AV  all  cracks  will  be 

subject  to  the  same  macroscopic  stress.  We  also  assume  that  the  material  is 

macroscopically  homogeneous,  so  that  a function  NCa^)  can  be  defined  as  the 

number  of  cracks  per  unit  volume  having  a critical  stress  equal  to  or  less 

than  a . The  critical  stress  of  a crack  is  defined  as  the  remote  stress  which 
c 

will  cause  fracture  when  applied  normal  to  the  crack  plane.  The  probability 

that  a crack  having  a critical  stress  in  the  range  a ^ to  CTc  + da^  exists  in 

volume  element  AV  is  then  AV  [dN(a  )/d a ] da  . 

c c c 

If  such  a crack  does  exist,  the  probability  that  it  will  fracture  depends 
on  its  orientation,  the  stress  state,  and  the  fracture  criterion.  We  assume 
there  is  a solid  angle  ft  such  that  fracture  of  a crack  will  occur  if  and  only 
if  its  normal  lies  within  ft.  This  means  that  if  the  normal  lies  within  ft, 
ag  > ac,  where  ag  is  the  effective  stress  corresponding  to  the  fracture 
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criterion  selected.  If  the  cracks  are  randomly  oriented,  the  probability 
that  a crack  will  fracture  under  the  applied  stress  E is  n(E,oc)/4u. 

Now  the  probability  of  failure  due  to  a crack  in  the  critical  stress 
range  dac  located  in  volume  element  AV  is  the  product  of  the  above  proba- 


bilities, i.e.. 


dN(ac) 


da 


c 


(2) 


Substituting  equation  (2)  into  equation  (1)  and  changing  sums  into  integrals 
we  obtain 


ps = exp  [’/dv/da c fr  If]  (3) 

which  is  the  result  given  previously  in  Reference  4. 

We  note  in  passing  that  NCa^)  is  independent  of  stress  state  and  depends 
only  on  the  material.  Since  a is  defined  as  the  stress  which  causes  fracture 
when  applied  normal  to  the  crack  plane,  N(ac)  can  be  converted  directly  into 
crack  size  distribution  when  KIc  is  known.  In  this  respect  it  differs  from 

9 10 

Weibull's  n and  the  g(S)  of  McClintock  et  al.  * These  functions  represent 
the  number  of  cracks  per  unit  volume  that  will  be  fractured  by  a particular 
applied  stress  state,  and  they  depend  on  the  stress  ratio.  More  precisely, 
they  are  essentially  equivalent  to  the  second  integral  in  Equation  (3). 

The  function  g(S)  is  the  same  as  N for  the  particular  case  of  hydrostatic 
tension,  because  then  is  4tt  or  zero  depending  on  the  stress  level. 

A consideration  of  for  simple  tension  serves  to  illustrate  the  depen- 
dence of  this  function  on  the  fracture  criterion  selected.  The  normal  stress 
criterion  assumes  an  effective  stress 
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a - a • o cos  0 
e n 


(4a) 


**  + *&]•*-*  f. 


oriented  crack  of  critical  stress  0 will  be  fractured  depends  on  both  the 


In  simple  tension  and  equibiaxial  tension  analytical  expressions  can 


be  found  for  ft  and  used  together  with  Equation  (3)  to  evaluate  P_-.  In  the 


general  case  we  find  ft  by  integrating  dft  over  the  range  in  which  0 > a 


One  way  of  accomplishing  this  is  to  integrate  over  the  entire  angular  range 


dvdftda  H(o  .0  ) 


1 when  a > a 


We  can  now  carry  out  the  integral  over  a first,  with  the  result 


In  the  case  of  the  normal  stress  criterion,  Equation  (4a),  0 is  substituted 


for  0 in  Equation  (8) . The  resulting  equation  is  completely  equivalent  to 


Weibull's  unproved  rule  for  treating  polyaxial  stress  states,  and  shows  that 


this  rule  is  valid  if  the  material  obeys  the  assumptions  listed  in  the  Batdorf 


Comparison  With  Exper iment 


The  theory  just  outlined  is  next  compared  with  the  results  of  two  experl 


ments.  The  first,  performed  by  Petrovlc  and  Mendiratta 


involved  the  crea' 


tion  of  semicircular  cracks  in  the  surface  of  hot  pressed  Sl.N, . By  finding 


out  how  the  average  fracture  stress  of  these  controlled  cracks  subjected  to 


four-point  bending  varies  with  crack  orientation,  a direct  test  of  the 


fracture  criteria  of  Equations  (A)  is  possible.  This  comparison  is  shown  in 

Figure  2.  The  theoretical  curves  are  obtained  by  solving  Equations  (A)  and 

(5)  for  o with  0£  taken  as  equal  to  the  measured  fracture  stress  when  0=0. 

It  should  be  noted  that  the  material  fractured  because  of  natural  cracks  at 

2 

stresses  above  approximately  500  MN/mm  so  that  only  the  data  for  0=0°, 
22.5°,  and  A5°  are  useful  for  our  purpose.  It  appears  that  the  second  most 
shear  sensitive  criterion,  Equation  (Ad),  agrees  best  with  the  data.  However, 
the  criteria  do  not  include  corrections  for  the  fact  that  these  were  surface 
cracks.  The  effect  of  a free  surface  is  to  enhance  the  influence  of  normal 
stress  on  fracture,  i.e.,  to  reduce  the  relative  importance  of  shear.  It  is 
possible  that  if  the  criteria  were  adjusted  to  take  this  into  account  Equa- 
tion (Ae)  might  turn  out  to  be  the  best. 

The  second  experiment  was  conducted  by  Giovan  and  Sines^  on  high 
strength  alumina.  In  this  experiment  the  statistics  of  failure  were  taken 
for  rectangular  plates  in  A-point  bending  and  for  disks  subjected  to  concen- 
tric ring  loading.  The  scatter  of  the  data  points  is  such  that  elaborate 
techniques  for  closely  fitting  the  data  appear  unwarranted.  Accordingly  we 
employ  Weibull's  two-parameter  representation  which  has  substantial  advantages 
in  computational  simplicity.  Among  other  things,  this  allows  us  to  treat 
bending  as  though  it  were  uniform  tension  applied  to  the  reduced  volume 
V*  = V/(m  + 1),  where  m is  the  Weibull  exponent. ^ 

We  assume  that  V’N(oc)  can  be  expressed  in  the  form 

V’N(o  ) = Kam  (9) 

e c 

Substituting  in  Equation  (3)  and  assuming  the  material  is  uniformly  stressed 
in  simple  tension  we  obtain 
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(10) 


Pg(o»0)  = exp 


[-  TS 


do  mKom  / 2ir  sin  9d0 
c c J 

o 


where  0g  is  the  crack  angle  for  which  = o^.  Integrating  we  obtain 


= exp 


Pg(a,0)  = exp  mK^*  (l-cos0e)  <£  ^dO£ j 


where 


ky  = mK  C (l-cos©^)  xm  ^dx 


(ID 


(12) 


and 


x = o /a 
c 


Similarly  for  equibiaxial  tension 


(13) 


r /*0  . 0 -» 

P(o,o)  = mK  J do^mKo™  J ecos  0d0j 


kBomJ  (14) 

where 

k = mK  f sin  0 xm  ^dx  (15) 

B Jo  e 

A convenient  way  to  compare  uniaxial  and  equibiaxial  fracture  statistics 
is  to  plot  &nPg  (0,0) /JlnPg(a,0)  = kg/ky  against  the  Ueibull  parameter  m.  This 
is  done  in  Figure  3.  It  is  evident  from  this  figure  that  Pg(o,0)  differs  most 
strongly  from  Pg(o,0)  for  the  case  of  shear-insensitive  cracks.  Equation  (4a). 
The  curve  for  this  case  is  the  same  as  that  shown  in  Figure  1 of  reference  8. 
The  other  four  curves  correspond  to  the  other  fracture  criteria  listed  in 
Equations  (4). 


These  curves  can  be  used  to  analyze  the  experimental  results  of  Giovan 


i 

t 
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and  Sines,  which  are  shown  in  Figure  4.  Either  of  the  curves  can  be  con- 
sidered as  given  and  used  to  calculate  the  other  theoretically.  Here  we  take 
the  equlbiaxial  data  as  given  and  by  conventional  methods  find  for  the  best- 
fitting Weibull  formula 

Pg  = exp  [-9.32  x 10-25O15’12] 

With  the  above  choice  of  m we  find  the  k ratios  for  the  various  fracture  cri- 
teria listed  in  Equations  (4).  The  corresponding  uniaxial  fracture  probabili- 
ties are  plotted  in  Figure  4.  It  is  evident  that  the  best  agreement  is 
obtained  using  Equation  (4e)  as  fracture  criterion. 

Discussion 

In  this  paper  two  equivalent  formulations.  Equations  (3)  and  (8),  are 
given  of  weakest  link  theory  in  which  the  links  are  considered  to  be  indi- 
vidual cracks.  The  fracture  criterion  and  the  stress  state  in  these  formula- 
tions are  arbitrary.  When  the  cracks  are  assumed  to  be  shear  insensitive, 
i.e.,  only  the  component  of  applied  stress  normal  to  the  crack  plane  contrib- 
utes to  failure.  Equation  (3)  reduces  to  the  theory  of  Reference  4 and  Equa- 

3 

tion  (8)  reduces  to  Weibull' s rule  for  polyaxlal  stress  states. 

A number  of  fracture  criteria  for  shear  sensitive  cracks  are  considered. 
These  criteria  are  formulated  in  Equations  (4b)  to  (4e).  Generally  speaking, 
penny-shaped  cracks  are  more  shear  sensitive  than  Griffith  cracks,  and  criteria 
resulting  from  strain  energy  release  rates  are  more  shear  sensitive  than  those 
based  on  the  assumption  that  fracture  occurs  when  the  local  stress  exceeds 
the  material  strength  at  some  point  on  the  surface  of  the  crack  cavity. 

The  theoretical  predictions  are  compared  with  two  sets  of  experimental 
data.  The  first  set  (Petrovic  and  Mendiratta)  tested  the  variation  in  fracture 

13 


250 


275 


300 


325  MPa 


Fie.  4.  Theory  Compared  with  Data  of  Giovan  and  Sines 


strength  of  surface  cracks  with  angle  of  inclination  to  the  tensile  axis. 

In  this  case  Equation  (4d)  appeared  to  be  in  best  agreement  with  experiment, 
and  Equation  (4e)  appeared  to  be  too  sensitive  to  shear.  However,  the  cri- 
terion were  not  adjusted  for  surface  effects  which  tend  to  reduce  shear  sen- 
sitivity, so  this  result  may  be  misleading.  In  the  second  set  of  experiments 
(Giovan  and  Sines)  Equation  (4e)  appeared  to  be  in  the  best  agreement.  Thus 
in  both  experiments  the  criteria  derived  from  strain  energy  release  rate  con- 
siderations worked  out  the  best. 

It  is  of  some  interest  that  in  the  second  set  of  experiments  a criterion 
exhibiting  greater  shear  sensitivity  than  any  considered  would  lead  to  still 
better  agreement  with  experiment.  The  evidence  for  greater  shear  sensitivity 
is  not  very  compelling  because  the  discrepancy  between  the  predicted  and 
measured  mean  fracture  stresses  in  simple  tension  is  well  below  the  standard 
derivation  that  would  be  calculated  for  this  mean  based  on  the  fact  that  there 
were  only  ten  test  specimens  for  each  stress  ratio.  However  theoretical  con- 
siderations also  suggest  greater  shear  sensitivity,  since  the  strain  energy 
release  rate  criteria  used  are  based  on  the  simplifying  assumption  that  cracks 
extend  in  their ‘own  plane,  an  assumption  leading  to  too  high  a fracture  stress. 
It  is  unfortunate  that  at  the  time  of  writing,  substantial  differences  of 

opinion  exist  concerning  how  to  calculate  the  stresses  that  cause  mixed  mode 
18—22 

fracture.  Until  it  becomes  clear  which  criterion  is  best  on  theoretical 

grounds,  experimental  data  cited  in  the  paper  suggest  it  may  be  satisfactory 

to  use  the  criterion  of  Equation  (4d).  It  has  the  advantage  of  simplicity, 

it  is  probably  conservative  when  uniaxial  data  are  used  as  the  basis,  and  it 

is  almost  certainly  preferable  to  the  normal  stress  criterion  employed  explic- 

4 

ltly  by  Batdorf  and  Crose  and  implicitly  by  Welbull.  In  any  event.  Equations 
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